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Abstract. - We show that disorder increases the persistent current of a half-filled one- 
dimensional Hubbard- Anderson ring at strong interaction. This unexpected effect results from 
a perturbative expansion starting from the strongly interacting Mott insulator ground state. 
The analytical result is confirmed and extended by numerical calculations. 



The interplay between disorder and electron-electron interactions is a challenging problem, 
at the heart of the understanding of recent important experimental findings. In particular, 
the measured persistent currents [1] in metallic mesoscopic rings are much larger than the 
theoretical prediction for noninteracting electrons [2], and the metallic phase appearing in 
two-dimensional electron gases at low-density [3, 4] is at odds with the scaling theory of 
localization [5] for noninteracting disordered systems. 

Given the complexity of treating disorder and interactions on the same footing, only lim- 
iting cases or simplified models have been solved. In such cases, the emerging conclusion is 
that both, disorder and a repulsive interaction suppress the mobility and persistent currents. 
This suppression is expected to be even more important in mesoscopic systems of reduced 
dimensionality, where the role of correlations is enhanced by the poor screening. 

These intuitive conclusions have been challenged in few-particle models, where it has been 
shown that a repulsive interaction, in the presence of strong disorder, may result in an en- 
hancement of the electron mobility [6,7]. Such delocalization effects are very dependent on 
the model and its particular dimensionality, and their connection with the experimentally 
relevant cases is still a matter of debate. 

In this work we show that disorder may also have the unexpected effect of favoring the 
zero-temperature persistent currents, in the case of strongly interacting electrons in one- 
dimensional half-filled Hubbard rings. This result is established from a perturbative expansion 
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starting from the strongly interacting limit, and numerical calculations based on the Density 
Matrix Rcnormalization Group algorithm [8,9]. Within a different system, namely a two- 
dimensional disordered Hubbard model, numerical indications of such a behavior for the finite 
temperature conductivity have recently been reported [10], and related with the disorder 
induced breaking of the particle-hole symmetry. 

A large variety of analytical and numerical tools have been developed and used for the 
study of interacting fermions in disordered one-dimensional systems. In particular, for spinless 
fermions, it has been numerically shown that repulsive interactions suppress the persistent 
currents, provided that the disorder is not too strong [11, 12]. Only for very strong disorder, 
an enhancement due to repulsive interactions occurs [13-15]. In this case, there is an important 
difference between the behavior of individual samples (exhibiting pronounced current peaks) 
and ensemble averages (showing a broad region of small enhancement for moderate values 
of the interaction). Once the spin is included in the description (Hubbard- Anderson model), 
renormalization group [16] and numerical approaches [17, 18] show that, outside half filling, 
repulsive interactions enhance the persistent currents in weakly disordered systems. Our work 
extends these results towards half filling and strong disorder, showing that the interaction 
induced increase of the persistent current is much more important than in the spinless case. 

We write the one-dimensional Hubbard- Anderson Hamiltonian as H = Hq + Hi with 

MM M 

Ho = U^ni^^ni,l + W^^v,ni,^ and = -i XI X! ("^Uci-i.tr + 4-i,<7Ci,<7) • (1) 

2—1 i—1 a i—1 a 

Hq accounts for interaction and disorder energy, and Hi for kinetic energy. The operator 
c^ ^ (c- ^) destroys (creates) a particle with spin a on site i, rii^a = ^c^ ^ is the occupation 
number operator, and the hopping element t = 1 sets the energy scale. The Vi are independent 
random variables, uniformly distributed within [—1/2; 1/2]. W and U denote disorder and 
interaction strengths, respectively. We concentrate on chains with an even number of sites M 
and N — M electrons (half filling). The ring topology results from the boundary condition 
Co = CMexpi$, where the phase $ = 27r(/)/$o accounts for a magnetic flux (/> threading the 
ring, and $o = hc/e is the flux quantum. 

The flux dependence of the many-body ground state energy E{^) is periodic with pe- 
riod 27r. In localized systems, the persistent current at zero temperature / = —dEo/d^ « 
{AE/2) sin $ can be characterized by the phase sensitivity AE = E{Q) — E{tt), measuring the 
difference between the ground-state energies at periodic ($ = 0) and anti-periodic ($ = tt) 
boundary conditions. We shall concentrate in the sequel on two related quantities: AE and 
the stiffness D = M| Ai<^|/2, which can be linked to the conductivity [19] and the conductance 
of the system [20]. 

At half fllling, strong repulsive interactions lead to a Mott- Hubbard transition and thereby 
to a strong suppression of the persistent current [13] . For f/ t, the kinetic energy part Hi 
represents a small perturbation to the dominating term Hq of the Hamiltonian (1), and it is 
possible to expand the ground-state energy in terms oit/U [15,21-23]. Since H is symmetric 
with respect to spin rotation and conserves the projection Sz of the total spin S onto the z- 
axis, we choose to work in the subspace which is characterized by 5*^ = 0, without restricting 
S. Within this subspace, the TV-body eigenstates \ipa) of Hq are given by products of on- 
site states, and completely specified by the functions i\} {a) that select the sites which are 
occupied by the A''/2 t(J,)-electrons {k = 1, 2, ....N/2) in the configuration a. They read 

(n/2 \ (n/2 \ 
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The corresponding many-body energies are Ea = + E^, with the disorder contributions 
E^ = WJ2k + ^i^a)) ' ^^'i interaction energies E^ = Uea- We note Ca the num- 

ber of doubly occupied sites. For U ^ W,t there is a gap separating the states without doubly 
occupied sites (subspace <S) from those having at least one doubly-occupied site (subspace V). 
The <S-states are Mott-insulator configurations (e^ = 0) and all have the same lowest en- 
ergy Eq' = W^^Vi, independent of their spin structure. This degeneracy is broken by the 
hopping part Hi of the Hamiltonian. Restricting ourselves to second order, we obtain an 
antiferromagnetic Heisenberg Hamiltonian as the effective model within the subspace S [24]. 
Whenever neighbouring sites (i and i — 1) have opposite spins, virtual hoppings i i — 1 —>■ i 
and i — 1— 1 lead to a fltix-independent coupling —2t'^/[U{l — {vi — Vi-i)'^W^ /U^)] 
between the two states that differ only in the orientation of the i and i — 1 spins. 

According to a theorem by Marshall, the non-degenerate ground state of the antiferromag- 
netic Heisenberg Hamiltonian (with energy Eq) can be written as |xo) = J2/3es //sIV'/s)) where 
the weights //3 can be chosen real and positive with < //? < 1 [24] . The effective model also 
contains diagonal matrix elements which are given by the sum of the above mentioned terms 
over all pairs of adjacent sites with opposite spin. The two configurations with alternating 
spins {i]}^^ = 2fc and i^'^' = 2fc — 1) have the lowest diagonal elements matrix elements and 
therefore the largest weights f^. 

We calculate the leading contribution to the phase sensitivity AE, using the difference of 
the higher order corrections for the ground-state energies at periodic (<!? = 0) and anti-periodic 
($ = tt) boundary conditions. The leading flux-dependent corrections to the ground-state 
energy appear in M^^ order [23] as 

p(M) ^ V /• f \Hl\lp-/i){^l''n\---\ V'7M - 1 ) (V'7M -i\Hl\^i3') 

Z^W ^Eo-E^,){Eo-E^,).-.{Eo-E^^_,) ' ^> 



71,72, ...jTm-i f3,0' 



/? and 0' run over all the <S-states and the 7's over the intermediate X>-states |V'7)- The 

numerator of the terms in (3) is composed of matrix elements (-^-y, |-ffi|V-'7(+i} of the hopping 
part of the Hamiltonian, thus non-zero contributions to E^^^ arise only when all subsequent 
virtual states can be connected by one-particle hopping processes [23]. The resulting 
correction to the ground-state energy can be expressed using sums over sequences A^^'^ ^ 
of one-particle hopping processes, starting at \tp(3'} and ending at \ipi3)- 

The dependence of the corrections (3) on the boundary condition enters via the hopping 
terms between sites 1 ^ M whose sign is reversed for anti-periodic boundary conditions. For 
periodic boundary conditions, all non-zero hopping elements are —t and the numerators in (3) 
are given by (— t)"signPf (A) signP|(A), with P|(|)(A) being the permutation of the positions 
for f (I) electrons on the ring, resulting from the sequence A. For anti-periodic boundary 
conditions, an additional sign (—1)'"' appears, where hh is the number of hops across the 
boundary 1 <-» M contained in A. 

Only sequences A with odd hy, yield energy corrections whose sign depends on the bound- 
ary condition. They are the only contributions to the phase sensitivity AE. The lowest order 
sequences having odd hh (and connecting one ground-state component l^p/s) from S to another 
one) are those which contain exactly one hop between each given pair of neighbouring sites. 
All hops must be in the same direction along the ring and the net effect is the transfer of 
one particle around the ring. This is why we work in order M, which gives the lowest order 
contributions with odd hy, = 1 [25]. Typical sequences contributing to AE^^^ are shown in 
Fig. 1. Diagonal contributions (/3 = /?') occur due to processes of type (a): only electrons 
with a given spin move while the electrons with the opposite spin orientation remain on their 
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Fig. 1 - (a,b): Schematic illustration of hopping sequences of order n = M contributing to the phase 
sensitivity (4), for the example M = 6. (c): The disorder dependence of |Ai5| for an individual sample 
with M = N = IQ. The data (symbols) are compared to the result of our perturbation theory (lines). 



initial site. Non-diagonal contributions {(3 ^ /?') can also appear. An example for which each 
of the electrons moves by one site is depicted in Fig. 1(b). Thus, the contributions to AE'-^^ 
are given by those hopping sequences which contain either M forward or M backward hops. 
Using (3), we can express 

as sums over all forward hopping sequences A^'f'^ ^ , with an additional factor of two accounting 
for the corresponding backward sequences. We have defined — Eq^ — Wd^^ , and extracted 
the dominant parametric dependence U^^^ from the denominator. All contributing sequences 
lead to a cyclic perturbation of the N/2 operators corresponding to electrons with a given 
spin direction in (2), and since the weights f/^ are all positive, the sign of the phase sensitivity 
at strong interaction is given by (—1)^^^, as in the non-interacting case [23]. 

In order to study the effect of disorder on the phase sensitivity at strong interaction 
{U ^ W), we expand (4) in powers of W/U. Up to second order, this yields 

^ 0,0' \ ^ ; ^T! ^ \i S Km J J 

(5) 

Since the matrix elements of the antiferromagnetic Heisenberg Hamiltonian are functions of 
{W/UY, the components fp of its ground state must be even functions of W/U. Up to second 
order, ^ » ff'' + ff^'W'^ /U"^ ■ The normalization condition E/3 = 1 implies ff^ ff'' = 
0, and since f^^"-' > wc must have positive and negative coefficients f^^p . The dominating term 
in (5), corresponding to the clean limit — > 0, exhibits an interaction- induced suppression of 
the persistent current oc U{t/U)^ . This suppression gets exponentially more pronounced as 
the system size increases. As in the spinless case [15], this can be interpreted as an interaction- 
dependent localization length ^{U) = l/ln{U/t). 

In contrast to the cases of lower filling or spinless fermions, where the linear correction in 
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W decreases the phase sensitivity [23], in the half- filled Hubbard- Anderson chain, the first- 
order correction in W/U vanishes exactly. We demonstrate this in the following way: a given 
sequence of hops can be characterized by the subsequent positions of doubly occupied sites 
and empty sites. Starting from this, one can construct a second sequence by exchanging the 
positions of the doubly occupied and empty sites. These two sequences have the same numbers 
of doubly occupied sites Ca, , and the coefficients dai have opposite sign. Their contributions 
to the first-order correction in W/U cancel each other, and the sum over all sequences of the 
linear term in the expansion (5) vanishes. The same argument applies to higher odd orders 
of the expansion, all of which vanish. 

Thus, the second order term of (5) determines the disorder dependence of the persistent 
current atU^ W. The first component of its prcfactor is a sum over M — 1 positive quantities 
of order one, and can be estimated to J^ii^ai/^ai)"^ ^ M per sequence. The second component 
is a sum over terms which are also of order one, but do not have a preferential sign. In the 
sum over all sequences, this yields, assuming random signs, a correction which scales with the 
square root of the number of sequences A^a ^ M\, and can thus be neglected as compared to 
the first component which scales oc A^a- The same holds for the second order term emerging 
from the expansion of the fp. Further, the fourth order correction (not explicitly written in 
Eq. 5) is dominated by a positive term yielding a prefactor of the order per sequence. 

This leads to the important conclusion that disorder increases the phase sensitivity (and 
therewith the stiffness) and the persistent current of strongly interacting electrons. Moreover, 
the relative importance of the corrections increases with the system size M. 

The increase is specific to half filling. At other fillings, the symmetry between doubly 
occupied and empty sites disappears, and the non-vanishing first order term decreases the 
phase sensitivity. This insight is consistent with the conclusion of Ref. [10] that particle-hole 
symmetry breaking leads to a disorder-induced increase of the low-temperature conductivity 
in two dimensions. 

In order to verify and extend the analytical result derived above, we performed extensive 
numerical calculations for the examples of A" = 10 particles (5 t;5 i) on M = 10 sites 
and A = 20 particles (10 t; 10 i) on M — 20 sites, using the Density Matrix Renormalization 
Group algorithm [8,9]. This algorithm allows to compute the many-body ground-state energies 
£^(0) and E{tt) for the interacting disordered system with sufficient precision to evaluate 
AE = E{0) - E{tt) [26]. For A' = 10 (20), we found negative (positive) AE for all values of 
U and W, consistent with the non-interacting case and the analytical result for very strong 
interaction oc (—1)^/^, in agreement with the result of Ref. [23]. 

The disorder dependence of \AE\ for AT = M = 10 particles is shown in Fig. 1 (c) for a 
typical individual sample (disorder realization) , at several large values of the interaction. The 
numerical results are very well described by Eq.(5) for small ratios of W/U. In order to have 
an agreement over a large range of disorder (like in Fig. 1 (c)), we need to go to fourth order 
in W. The parametric dependence \AE\ = {Aq/U'^){\^ A^iyV /Uf ^ A^(W /Uf ) (solid lines) 
yields the values Aq k, 390000, A\ w 2.33, and A^ « 9.36 (valid for all curves), consistent 
with the above presented estimations. 

In Fig. 2, we show the interaction dependence of the stiffness D for M = A = 10 (a) and 
M = A/' = 20 (b) for typical individual samples (lines) and ensemble averages (symbols), for 
different disorder values W^, and for the clean case W = {) (solid lino). The ensemble averages 
(of logZ?) are performed over 100 different samples (disorder realizations). At C/ = 0, the 
disorder leads to Anderson localization, and D is strongly suppressed by increasing disorder 
(note the logarithmic scale). While the interaction always reduc;c;s the stiffness in clean rings, 
consistent with a Luttinger liquid calculation [27] , a weak repulsive interaction U <t leads to 
an increase of the stiffness when disorder is present. Such an increase was predicted from a 
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Fig. 2 Interaction dependence of the stiffness for A4 — N — 10 (a) and M — N — 20 (b). Lines 
represent log D as a function of the interaction strength U, for different values of the disorder W, 
for an individual sample. Symbols in (a) correspond to the ensemble average (logD). The statistical 
errors are smaller than the symbol size. In the insets, strong-interaction data are plotted on log-log 
scale, confirming the power law D ~ f7-(^-i) (thick lines). 



renormalization group approach [16], away from half filling. Here, at half filling, the increase 
becomes very pronounced as the disorder increases. 

At large repulsive interaction U ^ W, in the Mott insulator limit, the behavior is radically 
different. The stiffness decreases strongly with the interaction, and, as shown in the insets of 
Fig. 2, the numerical results agree with the predicted power law oc 1/U^~^. While such an 
agreement is already obtained for individual samples, it persists for the average values. More- 
over, and in contrast to the Anderson insulator, the disorder increases the phase sensitivity of 
the Mott insulator. This enhancement can be very strong and become larger than an order of 
magnitude. At intermediate interaction, between the two previous limits {t U W), the 
phase sensitivity exhibits a maximum which becomes broader with increasing disorder and 
the disorder-induced decrease of its height is much weaker than for the Anderson insulator. 

Our analysis over two sizes confirms that the disorder-induced increase of the persistent 
currents becomes more pronounced (in relative terms) as the system size increases. In the case 
oi N = M = 10, where the ensemble average can be readily performed, we observe a similar 
behavior for individual samples and the average. Noticeable differences only appear at very 
strong disorder, but the analytically predicted enhancement of the phase sensitivity is robust 
with respect to impurity average. For the large system size M = A'' = 20, small peaks appear 
in the interaction-dependence of D, in individual samples at strong disorder. This structure 
is less pronounced than in the case of spinless fermions in strong disorder W = 9 [14, 15], 
where sharp peaks appear at realization-dependent values of U. It is also important to notice 
that the data of Fig. 2 exhibit, for all finite values of the disorder, a pronounced increase of 
the stiffness at intermediate U (with respect to the noninteracting case). This enhancement 
for the average can be larger than an order of magnitude, while for spinless fermions a small 
maximum appears only at very strong disorder [14,15]. 

In conclusion, we have investigated the effect of disorder on the persistent current of 
interacting electrons in one-dimensional half-filled chains, taking into account the spin degree 
of freedom. Within a systematic perturbative expansion of the phase sensitivity, we have 
unambiguously demonstrated that disorder increases the persistent current in the presence 
of strong interactions, in the Mott insulator limit. This is in striking contrast to the non- 
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interacting Anderson insulator limit, where disorder localizes the electrons and suppresses the 
persistent current. Such an effect can be related to the disorder-induced reduction of the 
energy gap between the Mott insulator ground state and excited states. Moreover, we have 
also shown that, for intermediate interactions between the two strongly insulating Anderson 
and Mott limits {t < U < W), the stiffness exhibits a pronounced max;imum, consistently 
with the results away from half filling. Such an effect might be a precursor of an even stronger 
tendency towards interaction-induced delocalization in disordered two-dimensional systems. 
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